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ABSTRACT 

Conaon r ock-f orning Minerals have theraal infrared spectral features that are 
Measured in the laboratory to infer coapos it ion . NASA non operates an airborne 
Daedalus scanner (TINS) that collects six channels of theraal infrared 
radiance data (8 to 12 ua) that nay be used to Measure these saae features for 
rock ident if icat ion . Previously, false-color coaposite pictures aade fron 
channels 1,3, and 5 and enittance spectra for saall areas on these iaages have 
been used to aake lithologic Baps. In this article, central wavelength, 
standard deviation and aaplitude of noraal curves regressed on the enittance 
spectra are related to coaposit ional infornation for crystalline igneous 
silicate rocks. As expected , the central wavelength varies systeaat ical ly 
with silica content and with nodal quartz content. Standard deviation is less 
sensitive to coapos it ional changes, but large values nay result fron nixed 
rock types within a pixel. Aaplitude varies with quartz content and also with 
adnixture of vegetation . Conpression of the six TINS channels to three iaage 
channels aade fron the regressed paraaeters nay be effective in iaproving 
geologic napping fron TINS data, and these synthetic iaages nay fora a basis 
for the reaote assessaent of rock coaposition . 


INTRODUCTION 


Common rock-forming minerals have thermal infrared spectral 
characteristics that can be used to infer rock compositions and 
construct lithologic maps from remotely sensed images. The 
strongest spectral features lie between 8 and 10 pun and are 
attributed to the fundamental vibrations of Si-0 Ee.g., Hunt and 
Salisbury, 1974; Hunt, 19S03. These reststrahlen bands are most 
pronounced in quartz. F'artial substitution of A1 , Mg or other 
elements for Si in a wide range of silicate minerals results in 
changes in the width, depth and wavelength of the band, as well 
as changes in its fine structure Ce.g., Lyon, 19653. For 
example, reststrahlen bands for quartz (SiO = ) are centered near 
8. 6-9. 2 aot; for albite (NaAlSia>Q 3 ) >-ajfc about 8. 7-9. 9 awi; and for 
olivine ( (Mg,Fe) 3 SiCU) at about 10.0-10.5 pun. Although rock 
spectra are not simple mixtures of mineral spectra, systematic 
differences in the emittance spectra of rocks do correspond to 
features in the spectra of the major constituent minerals. These 
differences are sufficient to produce distinct differences in the 
apparent temperatures measured remotely, at different wave- 
lengths. Thus multi spectral thermal measurements are a means of 
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assessing rock composition, and multichannel thermal infrared 
images are a means of extending this information regionally, to 
form a basis for lithologic maps Ce.g., Vincent et al . , 1972; 
Kahle and Rowan, 19801. 

Recently, NASA has acquired a Thermal Infrared Mul ti spectral 
Scanner (TIMS) that measures six channels of photon flux data in 
the spectral region 8-12 //m EPalluconi and Meeks, 19853. 
Previously, TIMS data have been displayed as false-color pictures 
made from three of the six image channels, enhanced to emphasize 
emittance information ESoha and Schwartz , 1978; Gillespie et al . , 
19863. This approach does not make use of information from the 
discarded image channels. In a different approach, thermal 
emittances may be calculated from each of the image channels and 
displayed as spectra [Kahle and Walker, 19843. Although in this 
approach no data are discarded, the obvious advantages of photo- 
interpretation cannot be exploited because the data are no longer 
displayed as pictures. If the emittance images themselves are 
displayed as false-color pictures, then three channels must still 
be discarded. 

This paper discusses a method that combines the virtues of both 
the pictorial and spectral analyses. I model the calculated 
emittance spectra with Gaussian functions, which have only three 
parameters: a central wavelength, a standard deviation, and an 
amplitude. In principle, these may be displayed using color 
pictures. Because the displayed parameters are calculated from 
all six measured values, more information is displayed than in 
false-color pictures made from the acquired data. In this paper , 
I attempt to show that significant lithologic information is 
contained in these calculated parameters, and that they vary 
systematically with conventional measures of composition. The 
three parameters thus form a basis for quantitative photo- 
interpretation and lithologic mapping. 


DATA ACQUISITION 


TIMS Data 

TIMS is a Daedalus airborne line— scanner operated by NASA on its 
NSTL Lear jet. TIMS acquires six channels of radiance (photon 
flux) data in spectral bands having peak sensitivities at 8.3, 
8.7, 9.1, 9.8, 10.4, and 11.3 /jm, respectively. Details of the 
instrument are available in Palluconi and Meeks C19853. 

During the summer of 1982 TIMS was flown over granitic rocks in 
the Sierra Nevada batholith, California, and again over the 
Josephine ophiolite in northern California and southern Oregon, 
acquiring data for a wide range of crystalline igneous silicate 
rocks. TIMS data acquired over Death Valley, California EKahle 
and Walker, 19843 were used to augment these data. TIMS was 
flown near noon to maximize the thermal radiance and hence the 
signal /noise ratio. Data were provided to JPL by NASA/NSTL as 
computer-compatible tapes. 
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Ancillary Data 

Thermal reflectance spectra were measured for representat i ve rock 
samples collected in the field, on freshly broken and natural 
weathered surfaces. This was done at JPL using an Anal ect 
Fourier Transform Infrared Spectrometer . Spectra were averaged 
for roughly 25 locations on each rock sample, depending on grain 
size, to ensure that they represented the bulk rock composition 
and could be related to TIMS data collected over outcrop areas of 
~400 m = . 

Estimates were made of the modal mineral abundances in selected 
samples by point counts of thin sections CGillespie and Abbott, 
19843. SiO = content was measured by electron microprobe analysis 
of glass beads fused from powdered rock samples Ee.g., Nash, 
1964; Baldridge, 19793. Results were checked against values 
calculated from the modal mineral abundances for the same 
samples. 


DATA PROCESSING 

TIMS Data 

Calibration and Correction... The TIMS images were corrected for 
systematic panoramic distortions using a nearest-neighbor algo- 
rithm (not interpolation) , to avoid further spectral mixing than 
was inherent in the data as acquired. The data were converted to 
units of photon flux, using the line-by-line internal blackbody 
reference or calibration values. Some additional aspects of 
required image processing are found in a companion paper 
CGillespie, 19863. After the data were corrected for instrument 
effects, emittance values were calculated for the six image 
channels, using a variation of the method of Kahle et al . C19803. 

Calculating Eaittances From TIMS Data... In order to calculate 
emittances from TIMS data, it is necessary to assume an emittance 
for one channel because there are seven unknowns (six emittances 
and the temperature) but only six radiance measurements. Based 
on this assumption, a model ground temperature can be calculated 
for each point in the image. Emittances for the other channels 
are then found from the measured photon flux ratioed to the flux 
predicted by Planck's Law for a blackbody at the model tempera- 
ture. In this way, emittance images for all six channels were 
constructed from the acquired radiance images. 

Kahle et al . C 19803 assumed the emittance in channel 6 was always 
0.93. This is valid for granites, but not for mafic and ultra- 
mafic rocks, which have low emittances near 11 yum. For these 
rocks the reststrahlen band lies at longer wavelengths than 
channel 1, and the emittance in channel 1 is more-or — less 
constant. Therefore, for granitic rocks the assumption of Kahle 
et al . C 19803 was used, but for mafic and ultramafic rocks the 
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emi ttance in channel 1 was taken to be 0.98, a value that was 
measured -for several samples on the Ana-lect spectrometer . In 
the future a more elegant solution would be to calculate the 
apparent blackbody temperature -for each TIMS channel; the 
channel -for which the maximum temperature was -found would be 
assigned an arbitrary emittance (e.g., 0.95), and the remaining 
emittances would be calculated as described above. It would be 
useful to create an image that reported the channel which had the 
maximum temperature, also. For this study, however, the 
simplified scheme in which either channel 1 or 6 was fixed was 
used . 

Atmospheric Effects ... The method of calculating emittances 
described above does not correct for atmospheric effects such as 
absorption and radiation by water vapor. While model corrections 
have been applied EKahle et al . , 19811, these depend on some 
knowledge of the air mass above the scene at the time of data 
acquisition. For the images used in this study, that vital 
information was lacking. 

However, recent investigations with JPL's Portable Field Emission 
Spectrometer have shown that vegetation is spectrally flat in the 
8-12 //m region. Any observed structure in the TIMS emittance 
spectra of vegetation is thus a measure of actual atmospheric 
effects. Therefore, atmospheric effects were reduced by normal- 
izing the calculated emittances to values found for vegetation in 
the same scene. This approach is valid for multiplicative mixing 
(e.g., translucent atmosphere) but not for additive mixing (e.g., 
radiating atmosphere) . It is valid only for a similar optical 
path through a single air mass; thus the target and the vegeta- 
tion must be at similar scan angles and elevations, and must be 
viewed through the same mass of air. Water vapor affects 
primarily TIMS channel 1, whereas atmospheric ozone affects 
channel 4. The effect on channel 1 appears to be stronger and 
more variable than for channel 4. 

Ultimately, it will be desirable to correct for both 
multiplicative and additive atmospheric effects. The iterative 
spectral mixing models of Adams et al . E19863 could probably be 
modified to do this. One of their advantages is that it is not 
necessary to completely describe the atmosphere in order to 
identify and remove its effects on the image spectra. 

Sampling the Image Data... For this study, six-valued emittance 
spectra were extracted from the processed images for sites 
identified on an interactive image analysis station at JPL's 
Image Processing Laboratory. Sites were selected based on field 
experience and existing geologic maps [Moore, 1963; Cater and 
Wells, 1953; Harper, 1980, 19843 to represent all major rock 
types encountered in the scenes studied. They were chosen 
because they contained outcrops of rock larger than the ~20-m 
nadir pixels, no vegetation, and little soil. 
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Laboratory Spectra 

Thermal infrared reflectance spectra were also obtained in the 
laboratory for rock samples collected at some sites identified in 
the TIMS images. K'ahle et al . [ 19843 have demonstrated that, for 
remote observations of many natural rock surfaces, emittance and 
reflectance are strongly anti correlated at thermal infrared wave- 
lengths, and are in fact related through Kirchoff's Law: 


£ ( = i - r(» 


( 1 ) 


where e<\) is the emissivity or emittance and r(\) is the reflec- 
tivity or reflectance at wavelength \ Ce.g., Lyon, 19653. Such 
is not the case for all samples, for example, powders and 
polished surfaces CConel , 19693. I calculated emittance spectra 
according to equation (1). TIMS emittance values were simulated 
from laboratory data by integrating the product of the calculated 
emittance spectra and the TIMS instrument response functions for 
the six channels. 


Data Regression 

All TIMS and simulated TIMS six— valued emittance spectra for the 
selected sites and rock types were fitted by Gaussian functions 
according to 


£ (\) 



( 2 ) 


where A is the amplitude, x= is the central wavelength and <r is 
the standard deviation. Note that in equation (2) the Gaussian 
function is actually related to r<\>, not e ( \) . The fitting was 
done by linear regression on the cumulative distribution function 
of (1-eCX)) values plotted, on a probability grid, against the 
peak wavelengths of the six spectral bands of TIMS. The central 
wavelength of the fitted function was found from the abscissa 
CX> value for which the regressed line had an ordinate value 
(probabi 1 i ty ) of 0.5. The standard deviation was found from the 
difference between and the wavelength at which the probability 
was 0.68. The third parameter, the amplitude A of the fitted 
function, cannot be found in the above manner. Instead, A is 
taken to be the maximum difference found among the six calculated 
emittance values. 
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RESULTS 


TIMS and Laboratory Spectra 

Figures 1 through 5 show laboratory reflectance spectra and TIMS 
emittance spectra -for -fresh and weathered silicate rocks. 
Laboratory spectra -for freshly broken samples are shown in Figure 
1, plotted with reflectance increasing downward, to simulate 
emittance. All the samples came from the study areas, except the 
Eureka Quartzite, which was from the Death Valley data CKahle and 
Walker, 19843. Reflectances range from IX to 79% with the maxi- 
mum value (minimum emittance) occurring at different wavelengths 
from sample to sample, from 8.6 to 10.3 yum. 

The orthoquartzite exhibited the greatest emittance (ref lectance) 
contrast , with a double reststrahlen band from 8 to 9 yum which is 
characteristic of Si0 = . This doublet is also prominant in the 
spectrum of the quartz-rich (28%) Bullfrog quartz monzonite, 
along with other bands attributable to feldspar and possibly clay 
minerals. These bands are found from 9 to 10 yum, and make the 
spectrum asymmetric. The maximum reflectance is * v 37%, much less 
than for the Eureka Quartzite. This is probably due to the lower 
quartz content of the Bullfrog pluton and to the abundance of 
feldspar, which has less intense reststrahlen bands than quartz . 

The spectrum of quartz— poor (7%) rocks, also from the Bullfrog 
pluton, shows the quartz bands subordinate to the feldspar bands, 
which have peak reflectances of only 12%. The spectrum of the 
hornblende diorite, which has only a trace of quartz , has even 
lower reflectances for "k < ' v 9 <um. The feldspar (plagioclase) 
band near 9.5 yum is still present, but the largest band is 
probably due to the hornblende, centered near ' v 10. 2 yum. 

Finally, the spectrum of the Josephine Peridotite has an intense 
band (38% contrast) near 9.5 yum , and a large shoulder near 11 jusn. 
The dominant band is attributable to serpentine minerals; the 
smaller features are probably due to olivine. 

Figure 2 contrasts laboratory reflectance spectra for freshly 
broken and natural weathered surfaces of the quartz— rich Bullfrog 
quartz monzonite and the serpentini zed Josephine Peridotite of 
Figure 1. Weathered surfaces of both rock types are less than 
half as reflective as the fresh surfaces in the reststrahlen 
bands. Surprisingly , the contrast in quartz bands is reduced 
proportionately more than the bands at longer wavelengths, for 
the quartz monzonite. We tentatively attribute this to the 
presence of a thin clay coat on the weathered surface. 
Reststrahlen bands for the clay minerals from 9 to 10 yum would 
keep contrast high there. 

Contrast for the weathered peridotite is reduced even more than 
for the quartz monzonite. A broad feature is still evident at 10 
yum, however. Apparently, the serpentine minerals are not 
preserved an the weathered surface, which consists primarily of 
Fe oxides, orthopyroxene, and minor olivine. 
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Figure 1. Laboratory reflectance 

spectra of fresh silicate rocks Figure 2. Laboratory reflectance 

(reflectance increases downward). spectra of fresh and weathered rocks 
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tance spectra calculated from TIMS for atmospheric effects. 
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Figure 2 emphasises that remotely sensed images show natural 
sur-faces that are not compositionally the same as the fresh 
bedrock. Nevertheless, weathering is not necessarily so intense 
that bedrock composition cannot be inferred. 

Figure 3 compares the TIMS emittance spectrum for an outcrop of 
the quartzose facies of the Bullfrog pluton to the laboratory 
reflectivity spectra of fresh and weathered surfaces. Clearly, 
the TIMS spectrum does not resemble the spectrum of the fresh 
surface, but it is similar to the spectrum of the weathered 
surface. Even so, the emittance for TIMS channel 1 appears to be 
somewhat too low. Perhaps this is due to incomplete correction 
for atmospheric water, which affects channel 1 most strongly. 
It is also possible that the disagreement arises from an 
incorrect choice of the model emittance for channel 6. To first 
order, this would introduce a gradient in the apparent spectrum, 
because of the non-linear nature of the Planck function. 

Figure 3 reinforces the findings of Kahle et al . C19843, and 
demonstrates that TIMS emittances and the laboratory reflectance 
spectra used in this study are in fact related through Kirchoff's 
law (equation 1). 

Figure 4 presents five TIMS emittance spectra for different 
silicate rock types. As in Figure 1, the Eureka Quartzite spec- 
trum has the highest contrast. The emittance for TIMS channel 1 
is lower than expected from the spectrum of the quartzite, which 
we attribute to atmospheric water. Kahle and Walker C 19843 did 
not correct the quartzite spectrum for atmospheric effects by the 
method described here, because of the absence of vegetation in 
Death Valley. The minimum emittance of 0.22 for the other chan- 
nels is considerably less than the value of 0.77 predicted by 
Kirchoff's Law for the fresh quartzite. Perhaps this is due to a 
thin film of desert varnish of the natural weathered surface. 
Even so, the quartzite spectrum has more contrast than those for 
the igneous silicate rocks. 

The spectra of the igneous silicate rocks show a decrease in 
contrast and an increase in the wavelength at which the minimum 
emittance is found as the composition becomes more mafic, as 
expected . 

Figure 5 shows TIMS emittance spectra for serpentinized Josephine 
Peri dot ite and for lateritic soils developed over the peri dot ite. 
The soil spectrum is nearly flat, a result corroborated by data 
from the JPL Portable Field Emission Spectrometer . The perido- 
tite spectra differ from the laterite and from each other primar- 
ily in channels 4 and 5 (9.8 and 10.3 ,um) . In Figure 5 it is 
evident that a false-color picture of the kind used in most early 
studies of TIMS data, based on channels 1, 3, and 5, would dis- 
criminate between the fresh and weathered rock poorly, because 
channel 4 was discarded. This emphasizes the need to examine all 
six image channels in lithologic mapping. 
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FRESH 6EOROCK (OUARRY) 

NATURAL BEDROCK EXPOSURE 

•••• LATERITIC SOIL OVER PERIOOTITE 


Figure 5. TIMS emittance spectra for Josephine Peridotite 
and lateritic soils in northern California. The emittance 
for Channel 6 was set to 0.95 for comparison to Fig. 4. 




Figure 6. Continuous laboratory reflectance spectra (dashed 
line), simulated TIMS emittance spectra (light solid line), 
and Gaussian functions fitted to TIMS spectra (heavy solid 
line). Left figure is for serpentinized peridotite; right 
figure is for quartz monzonite. Kol mogorov-Smi rnov test 
statistic is in the upper left-hand corner. Parenthetical 
values are calculated for reflectance vs V 1 . 
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Figure 5 shows that soils developed atop bedrock need not 
resemble the unweathered substrate at all. However, it should be 
possible to discriminate between soil and bedrock, and even to 
estimate the degree and weathering on bedrock outcrops. 


Fitting Gaussian Functions to TIMS Data 

Figures 6 through 9 describe the regression of normal functions 
to the emittance spectra and show the relationship among the 
three parameters of the fitted curves: A, and <r. Figure 6 

shows laboratory reflectance spectra, simulated TIMS data and 
the curves fitted to them for peridotite and quartz monzonite. 
For peridotite, the continuous spectrum is simple, and the 

assumption that it may be fitted by a Gaussian function appears 
to be realistic. The continuous spectrum for the quartz 

monzonite contains several distinct features, and will be less 
well described by a Gaussian function. However, the complexity 
is not evident in the TIMS spectrum, which is severely 

undersampled . 

The normality of the TIMS spectra was tested with the one-sample 
Kol mogor ov-Smi rnov test Ee.g., Bradley, 1968; Till, 19741, using 
the modified tables of Lilliefors £19671 that apply when the 

population is unknown. The test results show that the hypothesis 
that the TIMS spectra are normal cannot be rejected at the QO m A 
confidence level. Fitting curves according to wavenumber rather 
than wavelength did not appreciably improve the fit. The use of 
normal curves to describe the six-valued TIMS spectra is probably 
justified, although this might not be true for image data having 
higher spectral resolution and more channels. 

Figure 7 presents a variation diagram relating A (emittance 
contrast) and (central wavelength) for normal functions fitted 

to TIMS emittance spectra for various silicate rocks. The data 
loosely follow a trajectory in this diagram, according to 
lithology. A general dispersion of rock types with \ e= is 
more pronounced than with A. It appears that A decreases as 
increases for granitic rocks, but for mafic and ultramaf ic racks 
A is unchanged over a 0.6-;um range of These results are 

predictable from the laboratory data presented above. However, 
the scatter of data - especially along the A axis - may arise in 
part from natural variability of weathered surfaces in the TIMS 
data. In particular, minor amounts of vegetation, including 
lichens, will strongly affect A- 

Figure 8 shows the effect of spectral mixing between bedrock and 
vegetation (blackbody or graybody spectrum) in TIMS data. The 
emittance contrast in the spectrum for quartz monzonite is prog- 
ressively reduced with increasing vegetative cover. This effect 
is noticeable for as little as 5% cover and profoundly influences 
the parameter A. There appears to be little effect on the shape 
of the spectra; therefore, we expect cr and X= to be less sensi- 
tive than A to partial vegetation cover. 
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It is important to recognize that, in general , spectral mixing in 
thermal IR images can be quite complicated, even if it is limited 
to binary mixing between only two endmembers or components. Not 
only do the emittance spectra characteristic of the endmembers 
radiating to the sensor differ, but the temperatures at which 
they radiate may also. Because of the exponential temperature 
dependence of radiance in Planck's Law, the effect of temperature 
differences can be profound. At the least, it implies that a 
particular mixed spectrum can be achieved by quite a range of 
endmember concentrations or fractions. Further complications are 
introduced because the shape of the blackbody spectrum is also 
temperature-dependent . Thus it is non-trivial to estimate the 
fraction of vegetation cover, for example, from thermal infrared 
data. Despite this. Figure 8 demonstrates that it is possible to 
detect the presence of mixing in a straightforward manner . 

Figure 9 shows the relationship between the regression parameters 
«■ and Xc= for various silicate rocks. Although c ranges from 0.6 
to 0.9, there is no obvious correlation with \ t= . We note that 
is smallest for the monomineralic quartzite and for one of the 
peridotite samples. Granitic and especially mafic rocks have 
larger values of <r. Probably <r is larger for rocks containing 
several minerals, and hence for complicated spectra containing 
reststrahlen bands at two or more wavelengths. Thus <r, like A, 
may be sensitive to spectral mixing. 


Regression Parameters and Rock Composition 

Figures 10 through 12 show how the parameters A and X= vary with 
two conventional measures of rock composition: modal quartz and 
silica abundance. Figure 10 relates A and modal quartz 
abundance. As suggested from Figure 7, there is a pronounced 
positive correlation. The value of A ranges from ~0.07 for rocks 
with no quartz to ''0.16 for granitic rocks with 30X quartz. For 
the Eureka Quartzite, A is 0.22. The parameter A is clearly 
controlled in part by the intense quartz reststrahlen bands. As 
seen in Figure 8, in most image data A will also be controlled by 
vegetative cover, thus complicating its interpretation. 

Figure 11 illustrates the covariation of Xc and modal quartz. 
X= decreases from 9.7 to 8.5 jum as quartz content rises from 0 to 
100X. Half that variation occurs for quartz-free rock types. 
Ultramafic, mafic and granitic rocks and quartzite are all 
separable by x«=. The discrimination of the peridotites and 
quartz-free diorites shows that Xc is not sensitive to quartz 
alone. 

Figure 12 illustrates the correlation of X.= and silica abundance 
for the same rock types as Figure 11. It is important to empha- 
size that Xe is not responding to silica content per se; rather, 
it is responding to the succession of minerals that occurs with 
increasing silica in igneous rocks. Low-silica rocks contain 
olivine and pyroxene, both of which have reststrahlen bands near 
10 Aim. Rocks with more silica contain amphibole and feldspar. 
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Figure 11. Variation in (central 
wavelength) with modal quartz for 
TIMS emittance spectra of various 
silicate rocks. 
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Figure 12. Variation in (central 
wavelength) with silica content for 
simulated TIMS spectra of various 
silicate rocks, calculated from 
laboratory reflectance spectra. 


which have bands from 9 to 10 um- High-silica rocks have 
feldspar and quartz which has bands from 8 to 9 jum. Because 
silica content is a basis for the classif ication of igneous 
rocks, it is useful to note the indirect correlation of >«= and 
silica abundance. 


The simple relationship between X= and modal quartz or silica and 
the separability of common igneous rock types establishes 
regression of normal or other appropriate functions on TIMS 
emittance spectra as a useful tool for the construction of images 


for photointerpretation. 


particular 


strong inferences aboui 


rock type can be made from Xc so images of Xe will be an aid in 
lithologic mapping. Addition of A and «■ images, along with Xr , 
as elements in false-color pictures would increase the amount of 
information available to the photoi nterpreter . 


SUMMARY 

Six-channel emittance images were calculated from radiance images 
acquired in the spectral region 8-12 ,um by the NASA airborne 
Thermal Mul tispectral Infrared Scanner (TIMS) over the Sierra 
Nevada batholith and over the Josephine ophiolite in California. 
The apparent emittances were fitted by Gaussian functions 
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intended to model the reststrahlen bands in the emittance 
spectra. Thus the information contained in the original six- 
channel image was, in a sense, "compressed" to three channels 
representing the width, amplitude, and central wavelength of the 
Gaussian function. The width was greater for mul timineral ic 
rocks and appeared to be sensitive to mixing of mineral or rock- 
types. Contrast was especially sensitive to quartz content. The 
central wavelength was correlated to modal quartz abundances 
found from thin sections and to silica abundances measured by 
microprobe analyses of rock samples from selected outcrops 
representing a range of rock types. It was the best single 
measure of rock composition. Peridotites, hornblende diorites, 
granodiori tes and quartz monzonites were all separable by the 
regression parameters. 

Vegetative cover as low as 5X influenced the apparent 
emittances, and especially the contrast in the spectrum. Yet, 
because vegetation radiates like a blackbody or a graybody, 
spectral mixing of radiation from vegetation and a substrate does 
not generally obscure features of the emittance spectrum of the 
substrate, except for contrast. 

Thermal image data appear to be a valuable adjunct to field 
techniques in geologic mapping of silicate rocks. They augment 
aerial photographs and near-inf rared images in that they depict 
differences in silicate minerals composing the rocks, not just in 
the color index, the degree of oxidation, or the extent of 
hydration. 
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